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Abstract 

We report on a toy model for the glassy compaction dynamics of granular systems, 
introduced and investigated in collaboration with Anita Mehta and Peter Stadler. 
A stochastic dynamics is defined on a column of grains. Grains are anisotropic 
and possess a discrete orientational degree of freedom. Gravity induces long-range 
directional interactions down the column. The key control parameter of the model, e, 
is a representation of granular shape. Rational and irrational values of e correspond 
to very different kinds of behavior, both in statics (structure of ground states) and 
in low-temperature dynamics (retrieval of ground states). 
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1 Foreword 



In this contribution I have chosen to report on a series of works done in recent 
years in collaboration with Anita Mehta and Peter Stadler [1,2,3,4], where we 
have introduced and investigated in detail a columnar toy model for the com- 
paction dynamics of granular systems in the glassy regime. The main focus 
of this review is on the properties of this model which are of interest in the 
broader context of non-equilibrium dynamics and complex systems. We em- 
phasize the relationship between statics (structure of ground states) and low- 
temperature dynamics (retrieval of these ground states). Surprisingly enough, 
some of the key features of this model belong to Serge Aubry's favorites, such 
as quasiperiodicity and hierarchical melting. 
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2 Definition of the model 

A columnar model for the compaction dynamics of granular media in the glassy 
regime has been elaborated and investigated in a series of works [1,2,3,4]. The 
model consists in a column of N sites, each site being occupied by a grain. 
Grains are assumed to be anisotropic in shape. They take, for simplicity, two 
possible orientations, referred to as ordered and disordered. Sites are labeled by 
their depth n — 1, . . . ,N, measured from the top of the column. Orientation 
variables are defined by setting cr^ = +1 if grain number n is ordered, and 
an = —1 if grain number n is disordered. A configuration of the system is thus 
defined by the N orientation variables {cr„}. 

The collective dynamics of the orientation variables is chosen in order to model 
the compaction dynamics of granular systems under the influence of a small 
vibration intensity. The crucial effect of gravity is taken into account by con- 
sidering directional long-range interactions. Grain number n feels the effective 
weight of the whole piece of column above it (m = 1, . . . , n — 1). As a con- 
sequence, the upper grain (m) influences the dynamics of the lower one (n), 
whereas the lower grain has no influence on the dynamics of the upper one. 
In other words, causality acts both in space and in time. Hence there are no 
finite-size effects: a finite system made of grains behaves exactly as the N 
upper grains of a larger (finite or infinite) system. 

More specifically, we consider a Markovian stochastic dynamics defined by the 
following transition rates per unit time: 

w{an = -M ^ (7„ = -1) = exp 

< 

w{o-n = -1 ^ CTn = +1) = exp 



where T is a dimensionless measure of the vibration intensity, referred to 
as temperature, whereas the effective activation energy A„ and the effective 
ordering field /i„ represent the effect of all the grains above grain number n. 

In the most general model, A„ and /i„ depend on all the orientation vari- 
ables a^n for m = 1, ... ,n — 1. The simplest non-trivial model corresponds 
to [3,4] 

Xn^ An, hn^em~ -m:!;, (2.2) 
where m+ (resp. m~) is the number of ordered (resp. disordered) grains above 



Xn + hn 
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grain number n: 

2 n— 1 2 'n—l 



ml = ^Y.O^ + (Jk). m- = -^(l-<7fe), (2.3) 
^ fe=i ^ fc=i 



so that + m„ — n — 1. 

The shape parameter £ is the essential control parameter of the model. It has 
to obey £ > in order for the model to be non-trivial. It will be identified with 
the slope involved in the geometrical construction given below. Rational and 
irrational values of e, respectively corresponding to regular and irregular grain 
shapes, turn out to correspond to very different kinds of static and dynamical 
behavior. 

The amplitude A of the effective activation energy gives rise to a dynamical 
length 



T 
I' 



^dyn = (2.4) 



such that the Arrhenius factor of the rates at depth n reads 

Un = exp (^--f^ = exp (^-^j • (2-5) 



3 Ground-state structure 



As the dynamical rules of the model arc fully directional, they cannot obey 
detailed balance with respect to any Hamiltonian. 

There is, however, a well-defined concept of ground state. In the zero-tempe- 
rature limit we have indeed 

wjan = -1 ^ (7n = +1) _ / / oo if /l„ > 0, , . 

n:(an^+l^an^-l)'^^[T ) \0 if /i„ < 0. ^'"''^ 



We are therefore led to define a ground state of the system as a configuration 
where the orientation of every grain is aligned along its local field: 

, ( +1 a hn> 0, _s 
a„ = sign/^„ = |_^ if K<0 (3.2) 
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Table 1 

Recursive eonstruction of ground states by means of (3.2). 
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-1 
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-1 
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? 


? 



(provided hn 7^ 0). The condition (3.2) leaves the uppermost orientation cti 
unspecified, as hi vanishes identically. In the following, we assume for defi- 
niteness that the uppermost grain is ordered: 

ai = +1. (3.3) 



Ground states can be built in a recursive way (see Table 1). The number 
and the nature of ground states depend on whether e is irrational (only ai is 
unspecified) or rational (infinitely many orientations are unspecified). These 
two situations will be considered successively. 



3.1 Irrationals: Unique quasiperiodic ground state 



For irrational e, all the local fields hn are non-zero. Table 1 implies that they 
lie in the bounded interval 

-l<hn<e. (3.4) 



Let us introduce the following superspace formalism. Consider the integers 
(m~,m^) as the co-ordinates of points on a square lattice. We thus obtain a 
staircase-shaped line joining the points (mf , mf) = (0, 0), (m2 , r«2 ) = (0) 1); 
etc. Vertical steps correspond to ordered grains, whereas horizontal steps corre- 
spond to disordered grains. Equation (3.4) defines an oblique strip with slope e 
in the (m~,m"'") plane, which contains the entire line thus constructed (see 
Figure 1). 

A unique infinite ground-state configuration of grain orientations is thus gen- 
erated. This configuration is quasiperiodic. The above construction is indeed 
equivalent to the cut-and-project method of generating quasiperiodic tilings 
of the line. This approach, introduced by de Bruijn in the mathematical liter- 
ature [5] , became then very popular in the context of quasicrystals [6] . 
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Fig. 1. Geometrical construction of the quasiperiodic ground state of the model for 
the golden-mean slope (3.9). The two ways of going around the first cell, marked with 
a circle, correspond to the two possible choices for the orientation of the uppermost 
grain (after [4]). 

The ground state thus constructed can alternatively be described in analytical 
terms. We have indeed 

— n — 1 — m~ = 1 + Int((n — hn — g{nfl), (3.5) 



where the rotation number fl reads 

n = ^— (3.6) 

l + s ^ ^ 



and the hull function g{x) is a 1-periodic discontinuous sawtooth function: 



^ Frac(x — fl) 



In these formulas Int(a;), the integer part of a real number x, is the largest 
integer less than or equal to x, whereas Prac(x) — x — Int(x) is the fractional 
part of x (0 < Frac(a;) < 1). 

There are well-defined proportions of ordered and disordered grains in the 
ground state: 

/+ = Q = ^_, /_ = i_Q=^_. (3.8) 
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This geometrical construction is illustrated in Figure 1 for the most familiar 
of all irrational numbers [7], the inverse golden mean: 

£ = $-1 = 1 = 2-$=^, ^ ^ ^ + 1 ^ 1.618033. (3.9) 

The corresponding configuration is given by a Fibonacci sequence [5,6,7]: 

Wn} = H \ \ \ \ \ \ \ \ 



3.2 Rational e: Degenerate ground states 



For a rational slope 

£ = -, i.e., Q = ^— , (3.10) 
q p + q 



in irreducible forms {p and q are mutual primes), some of the local fields /i„ 
generated by Table 1 vanish. The corresponding grain orientations (j„ remain 
unspecified. 

This feature of rational slopes is clearly visible on the geometrical construction. 
Figure 2, corresponding toe = 2/3, shows that some of the lattice cells, marked 
with circles, are entirely contained in the closed strip (3.4). Consider one such 
cell. The broken line enters the cell at its lower left corner and exits the cell 
at its upper right corner. It can go either counterclockwise, via the lower right 
corner, giving cr^+i = —1, (7„+2 = +1, or clockwise, via the upper left corner, 
giving (Jn+i = +1, c„+2 = ~1- Each marked cell thus generates a binary 
choice in the construction. This orientational indeterminacy occurs at points 
such that n is a multiple of the period p + q, equal to the denominator of the 
rotation number Q,. 

The model therefore has a non-zero ground-state entropy (zero-temperature 
configurational entropy, complexity [8]) 

E = J^ (3.11) 

p + q 



per grain. Each ground state is a random sequence of two well-defined patterns 
of length p + q, each of them made of p ordered and q disordered ones, so 
that (3.8) still holds for each of the ground states. The patterns only differ by 
their first two orientations. The first cases are listed in Table 2. 
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Fig. 2. Geometrical construction of the ground states of the model for the rational 
slope £ = 2/3. The marked cells, entirely contained in the strip, are responsible for 
the non-zero configurational entropy (after [4]). 

Table 2 

Patterns building up the random ground states for the first rational values of e. 
The second example with period p + q = 5 is illustrated in Figure 2 (after [4]). 
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4 Zero-temperature dynamics. Retrieval of ground states 

We have already anticipated in the previous section that the stochastic dy- 
namics of the model simplifies in the zero-temperature limit. Equation (3.1) 
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indeed becomes the deterministic rule 



a. 



n 



sign /i, 



'n 



(4.1) 



(again provided /i„ 0). 

The dynamical length ^dyn introduced in (2.4) is assumed to have a non- 
trivial zero-temperature limit. Most of the interesting features of the model 
are already present in the limiting case where ^dyn = 00 ■ 



J^A Irrational e, infinite ^dyn-' Ballistic coarsening 

For irrational £, the rule (4.1) is always well-defined, as the local fields /i„ never 
vanish. We assume that the system is initially in a disordered state, where each 
grain is oriented at random: (T„ = ± with equal probabilities, except for the 
uppermost one, which is fixed according to (3.3). 

Consider first the limiting situation where ^dyn is infinite. The zero-tempe- 
rature dynamics is observed to drive the system to its quasiperiodic ground 
state. This ordering propagates down the system from its top, via ballistic 
coarsening. At time t, the grain orientations have converged to their ground- 
state values, given by the above geometrical construction, in an upper layer 
whose depth is observed to grow linearly with time: 



whereas the rest of the system is still disordered. This phenomenon is similar 
to phase ordering, as order propagates over a macroscopic length L{t) which 
grows forever. It is however different from usual coarsening, as the depth of 
the ordered region grows ballistically, with a well-defined e-dependent ordering 
velocity V{e), instead of diffusively, or even more slowly [9]. 

The ordering velocity obeys the symmetry property V{e) — V{l/e). Besides 
this property, its ^-dependence can only be investigated numerically. The or- 
dering velocity is observed to vary smoothly with e (although it is only defined 
for irrational e), and to diverge as V{e) ~ 1/e as £ — > 0. Figure 3 shows a plot 
of the inverse ordering velocity against £, for < £ < 1. 



L{t) ^ Vt, 
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e 

Fig. 3. Plot of the inverse ordering velocity 1/V of zero-temperature coarsening 
dynamics at infinite ^dyn, against the irrational slope e, for < e < 1. 

4-2 Irrational e, finite ^dyn-' Crossover to logarithmic coarsening 



In the general situation where ^dyn is finite, but large at the microscopic scale 
of a grain, the ballistic coarsening law (4.2) is modified in a minimal way, 
by taking into account the dependence of the rate (2.5) on depth: dL/dt ~ 
V exp(— L/^dyn)- We thus obtain 



L{t) ^ ^dyn In 1 + 



Cdyn, 



(4.3) 



This result exhibits a crossover between the ballistic law (4.2) for Vt -C idyn 
and the logarithmic coarsening law 



L{t) ^ ^dyn Int 



(4.4) 



for Vt > ^dyn- 

Equation (4.3) has been checked against the results of numerical simulations, 
for the golden-mcan slope. Figure 4 shows a scaling plot of L{t) corresponding 
to ^dyn = 50 and 100, together with the prediction (4.3), with no adjustable 
parameter. The ordering velocity V ^ 2.5S is taken from the data of Figure 3. 
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dyn 



Fig. 4. Scaling plot of L{t) /^dyn against t/^dyn for zero-temperature coarsening 
dynamics with the golden-mean slope. Symbols: numerical data. Full line: predic- 
tion (4.3), with V = 2.58. 

4-3 Rationale: Fluctuating steady state with anomalous roughening 

We now turn to zero-temperature dynamics for rational e. The updating 
rule (4.1) is not always well-defined, as the local fields /i„ may now vanish. In 
such a circumstance, it is natural to choose the corresponding orientation at 
random: 



The zero-temperature dynamics defined in this way therefore keeps a stochas- 
tic component. 

For concreteness, let us focus our attention onto the limiting situation where 
^dyn is infinite, and onto the simplest rational case, i.e., e — 1. Equation (2.2) 
for the local fields reads 



a. 



n 



' +1 if hn>0, 

< ±1 with prob. 1/2 if /i„ = 0, 
-1 if hn<0. 



(4.5) 



n-l 




(4.6) 



m=l 
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Fig. 5. Log-log plot of the local field variance = {h'^) against depth n, for 
zero-temperature dynamics with e = 1. Fnll line: numerical data. Dashed line: fit 
to asymptotic behavior, leading to (4.7) with A ~ 0.83 (after [3]). 

It turns out that the zero-temperature dynamics is not able to retrieve any of 
the degenerate ground states. The system rather reaches a non-trivial fluctuat- 
ing non- equilibrium steady state. This steady state is unique, i.e., independent 
of the initial configuration of the system. It is reached after a microscopic 
time. 

This steady state is characterized by unbounded, albeit subextensive fiuctu- 
ations of the local fields. Figure 5 indeed demonstrates that the local field 
variance grows as 



The exponent 2/3 of this anomalous roughening law can be predicted by the 
following line of reasoning. Viewing the depth n as time, and the local field hn 
as the position of a fictitious random walker, the noise in this walk originates 
in the earlier times m such that hm = 0. The probability for to vanish is 
of the order of 1/Wm- The effective diffusion coefficient at time n therefore 



scales as ~ {^/n)J2m=i ^ hence ~ nD^ ~ n/Wn, i.e.. 



Another way of characterizing this fluctuating steady state consists in look- 
ing at the probabilities p(C) of all the microscopic conflgurations C. Figure 6 
shows a plot of the normahzed probabilities 2^^p(C), against the 2^^ = 4096 



W^ = {hl)^An'/\ 



(4.7) 
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Fig. 6. Plot of the normalized probabilities 2^^p(C) of the configurations of a 
column of 12 grains in the zero-temperature steady state with £ = 1, against the 
configurations C in lexicographical order. Open circles: degenerate ground states 
(after [4]). 

configurations C of a column of 12 grains for £ = 1, sorted according to lex- 
icographical order (read down the column). This plot exhibits a rugged and 
seemingly fractal structure with self-similarity at all scales. The most fre- 
quently visited configurations are the 2^ = 64 degenerate ground states of the 
system, shown as open circles. 

The knowledge of all the steady-state probabilities p{C) gives access to the 
entropy of the fluctuating steady state, defined by means of the Boltzmann 
formula 

5 = -^p(C)lnp(C). (4.8) 
c 

An estimate for the entropy S can be derived by using the main feature of the 
zero-temperature steady state, i.e., the roughening law (4.7). Think again of 
the depth n as a fictitious discrete time, and of the local field hn as the position 
of a random walker at time n. For a free lattice random walk of n steps, one has 
{h^) = n, and the entropy reads Sfiat = nln2, as all configurations are equally 
probable. In the present situation, the entropy S is reduced with respect to 
S'flat, because (hi) = <^ n. The entropy reduction can be estimated as 
= 5'flat — S ^ J2m=i n^^^ ■ Figure 7 shows a plot of the entropy 

reduction AS against n. A reasonable agreement with the above estimate is 
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Ql ^ ^ . ^ . ^ . 1 

5 10 15 20 

n 

Fig. 7. Plot of the entropy reduction AS" = S'fjat — 5" in the zero-temperature 
steady state with e = 1, against depth n. Symbols: numerical data. Full line: fit 
AS" = (621nn + 53)10"^ ra^/^ (after [4]). 

found. 

The total entropy S of the fluctuating steady state is therefore equal to S'flat = 
nln2, up to a small subextensive reduction A^" of order n}/^. The equahty 
of the extensive part of S with S^at can be viewed as a manifestation of the 
flatness hypothesis of Edwards [10]. 



5 Low-temperature dynamics 

We now turn to the investigation of the dynamics of the model for a low but 
non-zero temperature T. We consider for simplicity the case of an infinite ^dyn- 

If the slope e is rational, the rule (4.5) is already stochastic at zero tempera- 
ture, so that no qualitatively novel effect appears at low temperature. 

We therefore focus our attention onto the case of an irrational slope e. We recall 
that the zero-temperature dynamics (4.1) drives the system to its quasiperiodic 
ground state, where each orientation is aligned with its local field, according 
to (3.2). For a low but non-zero temperature T, there will be mistakes, i.e., 
orientations (T„ = — sign/i„ not aligned with their local field. Equation (3.1) 
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suggests that the probabihty of observing a mistake at site n scales as 

n(n) ^ exp (-^] . (5.1) 



Hence the most fragile sites n, such that the local field hn is relatively small 
in the ground state ~ T -C 1), will be preferred nucleation sites for 

mistakes, and thus govern the low-temperature dynamics. 

These nucleation sites can be located as follows. Equation (3.5) shows that 
the local field /i„ is small when nQ, is close to an integer m. The latter turns 
out to be m = m^. Indeed 

r 

nVL = m + 5^K = pr, {VL-l<5 <n). (5.2) 



The most active nucleation sites are therefore in correspondence with rational 
numbers m/n which are the closest to the irrational rotation number Vt. Find- 
ing these rational approximants is a well-defined problem of Number Theory, 
referred to as Diophantine approximation [7]. 

Let us illustrate this on the example of the golden-mean slope (3.9). We are led 
to introduce the Fibonacci numbers [6,7], defined by the recursion formula 



+ Ffe_2 (Fo = 0, Fi = l). (5.3) 



We have alternatively 



(6.4) 



The leading nucleation sites are the Fibonacci sites n — F^. We have m — 
m+ = Fk-2, m- = Ffe_i, and = (-)^$-(^-i), so that 

/ 2$ \ 

n. = n«).exp - , (5.5) 



We can therefore draw the following picture of low-temperature dynamics. 
Mistakes are nucleated at the Fibonacci sites F^, according to a Poisson pro- 
cess, with exponentially small rates proportional to the 11^. They are then 
advected with constant velocity V ~ 2.58, just as in the zero-temperature 
case. The system is ordered according to its quasiperiodic ground state in an 
upper layer (n < A/'(t)), while the rest is disordered, somehow like the zero- 
temperature steady state for a rational slope. The depth Af{t) of the ordered 
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Fig. 8. Plot of the instantaneous depth J\f{t) of the ordered layer against time 
t, for the golden-mean slope, and temperature T = 0.003. Dashed lines: leading 
nucleation sites given by consecutive Fibonacci numbers (bottom to top: Fu = 89, 
Fi2 = 144, Fi3 = 233) (after [3]). 

layer, given by the position of the uppermost mistake, is therefore a natural 
collective co-ordinate describing low-temperature dynamics. 

Figure 8 shows a typical sawtooth plot of the instantaneous depth A/'(t), for a 
temperature T = 0.003. The leading nucleation sites are observed to be given 
by Fibonacci numbers. 

The system thus reaches a steady state, characterized by a finite ordering 
length (jV), which diverges at low temperature, as mistakes become more and 
more rare. The law of this divergence can be predicted by the following argu- 
ment. The most active nucleation Fibonacci site is such that the nucleation 
time l/IIfc is comparable to the advection time to the next nucleation site 
Fk+i, i.e., (Ffc+i —Fk)/V Fk/{^V), hence UkFk ~ ^V. Less deep sites have 
too small nucleation rates, while the mistakes nucleated at deeper sites have 
little chance to be the uppermost ones. Equations (5.5) then yields 



2$ " 



1, 



(5.6) 



hence 



{^f) ^ ^ ^ (5.7) 
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The estimate (5.6) correctly predicts that F12 = 144 is the most active nucle- 
ation site at temperature T = 0.003 (see Figure 8). 

Both the existence of a well-defined sequence of preferred nucleation sites, 
and the divergence of the depth of the most active nucleation site as (jV) ~ 
|lnT|/T, are very reminiscent of the phenomenon of hierarchical melting, put 
forward by Serge Aubry and collaborators [11]. This effect manifests itself 
as peaks in the specific heat at low temperature in some incommensurate 
modulated solids, each peak being in correspondence with a temperature- 
dependent spatial scale, where the system is more fragile and preferentially 
melts. 



6 Outline 

In a series of recent works [1,2,3,4] we have introduced and investigated in 
detail a columnar toy model for granular compaction. The minimal model 
of interest described in this report has one essential parameter, the shape 
parameter e, which gives the slope of the geometrical construction of ground 
states within the super space formalism. 

Rational and irrational values of e correspond to very different kinds of be- 
havior, both in statics (structure of ground states) and in low-temperature 
dynamics (retrieval of these ground states). 

For irrational values of e: 

• The model has a unique quasiperiodic ground state. 

• Zero-temperature dynamics leads to a fast retrieval of this ground state by 
means of the ballistic growth of the thickness of the ordered top layer. 

• A crossover to a logarithmic growth takes place when this depth becomes 
comparable to the dynamical length. 

• Low-temperature dynamics involves the nucleation of mistakes at preferred 
nucleation sites. This phenomenon is analogous to hierarchical melting. 

For rational values e — p/q: 

• The model has extensively degenerate ground states and therefore a positive 
configurational entropy. 

• Ground states consist of random sequences of two well-defined patterns 
consisting of p + g grains. 

• Zero-temperature dynamics is not able to retrieve any of the ground states. 

• The system is rather driven into a non-trivial fluctuating steady state char- 
acterized by unbounded, albeit subextensive fluctuations in the local fields. 
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